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Elastic scattering observables (differential cross section and analyzing power) are calculated for 
the reaction ^He(p,p)^He at projectile energies starting at 71 MeV/nucleon. The optical potential 
needed to describe the reaction is derived by describing ^He in terms of a *He-core and two neutrons. 
The Watson first order multiple scattering ansatz is extended to accommodate the internal dynamics 
of a composite cluster model for the ®He nucleus scattering from a nucleon projectile. The calcu- 
lations are compared with the recent experiments at the projectile energy of 71 MeV/nucleon. In 
addition, differential cross sections and analyzing powers are calculated at selected higher energies. 
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I. INTRODUCTION 

Traditionally, differential cross sections and spin observables played an important role in either determining the 
parameters in phenomenological optical models for proton-nucleus (p-A) scattering or in testing the accuracy and 
validity of microscopic models thereof. Specifically, elastic scattering of protons and neutrons from stable nuclei lead 
to a large body of work on optical potentials in which the two-nucleon interaction and the density matrix of the 
nucleus were taken as input to ab-initio calculations of first order optical potentials, either in a KMT or Watson 
expansion of the multiple scattering series, for which the primary goal was a deeper understanding of the reaction 
mechanism. 

For exotic nuclei the theoretical emphasis is somewhat shifted. A goal of microscopic reaction theory is a construction 
of the scattering observables based on well-defined dynamical and structure quantities in order to, for example, examine 
structural sensitivities. Investigating the structure of halo nuclei, specifically ^He, has already inspired a large body 
of work including few-body models l]. Green's Function Monte Carlo (GFMC) methods [2|, and no-core shell model 
calculations [3], so that ground state properties of ^He appear to be quite well understood. 

Recently, elastic scattering of ^He off a polarized proton target has been measured for the first time at an energy 
of 71 MeV/nucleon [3, la]- The experiment finds that the analyzing power becomes negative around 50° which is not 
predicted by simple folding models for the optical potentials [6, 7]. The same calculations nevertheless describe the 
differential cross section at this energy reasonably well. This apparent "Ay problem" conveys the inadequacy of using 
the same models which describe p-A scattering from stable nuclei for reactions involving halo nuclei. The obvious 
difference is the nuclear structure. While the typical stable nuclei for which folding models are very successful are 
mostly spherical, ^He can be understood in few-body models as a three-body system consisting of 2 neutrons (n) and 
a *He core. Implementing this three-body structure in a cluster model, specifically in a reaction calculation for proton 
(p) scattering off ^He, was pioneered in Ref. Q for calculating the reaction p-l-^He at 717 MeV/nucleon, an energy at 
which the authors could employ the Born approximation. Based on the KMT formulation for the optical potential 
and more realistic wave functions for ^He, differential cross sections and analyzing powers were calculated in Ref. [9| 
at 297 MeV/per nucleon. For describing the differential cross section and the analyzing power at 71 MeV/nucleon, 
Ref. 4] suggested a "cluster-folding" calculation having an explicit a-core described by a phenomenological p+'^He 
optical potential. 

In this work we want to extend the development of Ref. [8| by incorporating the cluster structure in an optical 
potential for the reaction p-|-^He in such a way that the transition amplitude can be iterated to all orders (non-Born 
approximation). Our derivation of the optical potential is based on the Watson formulation for the multiple scattering 
theory, which not only allows to treat proton and neutron contributions to the structure separately (1^ , but also lends 
itself naturally to taking into account the contributions of the a-core and the two neutrons. The construction of an 
optical potential in which the separate contributions from the clusters are treated in a consistent fashion is achieved. 

This article is organized as follows. In Section II we first present a short summary of the Watson optical potential for 
stable nuclei, and then extend this paradigm to the ^He nucleus consisting of an a-core and two neutrons. In Section 
III we present our calculation for ^He-|-p at 71 MeV/nucleon as well as at several higher energies and discuss their 
implications. Our conclusions are presented in Section IV. Three Appendices are devoted to the explicit derivation of 
the first order optical potential, the transformations between the different coordinate systems used in our calculations, 
and the calculation of the correlation densities between the clusters. 



II. THE FOLDING CLUSTER MODEL 

In order to derive a cluster ansatz for the target (projectile), and show how it can consistently be incorporated into 
a folding approach for the optical potential, we will for the convenience of the reader give a summary of the essential 
ingredients and underlying assumptions. 



A. The Watson Optical Potential for Single Scattering 

Let H = Hq + V he the Hamiltonian for the nucleon- nucleus system (A-l-1 body system), where the interaction 
V — J2i ^Oj consists of all two-nucleon interactions voi between the projectile ("0") and a target nucleon ("i"). The 
free Hamiltonian is given by Hq = h^ + Ha, where ho describes the kinetic energy of the projectile, while the target 
Hamiltonian Ha satisfies Ha\'^a) — Ea\'^a), with |^a) being the ground state of the target. 

The transition amplitude for the scattering of the projectile from the target is then given by a Lippmann-Schwinger 



equation, T — V + VGqT, where the propagator Go is an (A+1) body operator given by 

Go{E)^{E~ho^HA + ie)-\ (1) 

with E being the total energy of the system. One way to tackle the many-body scattering problem is the spectator 
expansion 11], which writes the transition amplitude as T = X]i=i '^oi, so that 

To^^Vo^+Vo^Go{E)T. (2) 

This allows the sum of all interactions between projectile '0' and nucleon 'i' by a formal reordering of the multiple 
scattering series according to Watson, 

Toi = tqi + ToiGo{E)^ To J , (3) 

where 

Toi = vqi + vqi Go {E)Tot . (4) 

The term foi of Eq. Q only involves the interaction between pairs, namely particles '0' and 'i', whereas the propagator 
Go{E) is still an (A+l)-body operator. The multiple scattering series of Eq. Q can directly serve as starting point 
for constructing the transition amplitude for elastic scattering as shown in Refs. [8, 9, 12]. 

When focusing on elastic scattering, the projection P onto the ground state \^a) is introduced so that [Go{E),P] — 
0. Here we define P — /^ \ \j, A and P+Q — 1, where Q projects onto the orthogonal space. This allows the separation 
of the transition amplitude into two pieces, 

T = U + UGoiE)PT 

U = V + VGo{E)QU, (5) 

with U being the optical potential operator. The transition operator for elastic scattering may then be defined as 
T^i = PTP, so that 

T,i = PUP + PUPGo{E)PT,i (6) 

is a one-body integral equation. Of course, it requires the knowledge of PUP, which has to contain the complete 
information about the many-body character of the problem. The formulations for the transition matrix for elastic 
scattering, given in Eqs. ([3]) and ([5]), are equivalent though truncations in the expansions are not. 
The first order term of U can be defined as 



[/ = ^C/,«^ro., (7) 



with 

Toi = voi + voiGo{E)QToi. (8) 

Because of the appearance of the projection operator Q in Eq. ([7]), the quantity Tot can not yet be related to a two 
nucleon interaction. Defining a transition operator foi, according to Eq. (|4]), allows the explicit relation to Toi'. 

Toi = vot + voiGo{E) [P + Q] foi 

= To^ + TorGo{E)Pfo., (9) 

so that one obtains the exact relations [Ifll 

Toi = Toi - ToiGo{E)PTot = Toi - ToiGo[E)PToi. (10) 

Taking into account the iso-spin character of the target nucleons instead of just summing over A nucleons can be easily 
done by sphtting Eq. ([7]) into two parts under the assumption that the projectile '0' is a proton, Eq. d?]) becomes: 

Up -jz<f + E-oT ^ u^ + t/p^ (11) 



where the integral Eq. (TTO]) has to be solved separately for proton-proton (pp) and neutron-proton (np) interactions. 
This clearly indicates that the optical potential for the scattering of a proton {Up) from a target nucleus differs from 
the optical potential for the scattering of a neutron ([/„) from the same target. Moreover, and more important for 
the present considerations, the folding with the proton and neutron density matrices is cleanly separated, which is 
not the case if one uses an optical potential in the formulation of Kerman-McManus-Thaler (KMT) ISJ. A numerical 
study for p +^^ Li between a KMT formulation and Watson expansion of Eq. ([3]) has been carried out in Ref. [l2] 
with the finding that truncations at the same order of the series though being similar at small momentum transfer 
show differences at the higher momentum transfers. This should not be surprising when considering the nonlinear 
relation of the free two-nucleon t- matrix, Eq. (|10[) . to the quantity roi entering the optical potential which is the 
driving term of the final scattering integral equation, Eq. ©. For the n-\- d system, for which exact solutions of the 
Faddeev equations exist, optical potentials for elastic scattering were constructed in |14l | indicating that first order 
approximations are only valid for smaller momentum transfers. 

The propagator of Eq. ((!]) needs to be examined in more detail since it still is an (A-l-l)-body operator. Only if 
the target Hamiltonian Ha is approximated by a c- number, Gq{E) becomes a one-body operator and Eq. ([4]) can be 
identified with the free nucleon-nucleon (NN) t-matrix at an appropriate energy E. Reducing Ha to a c-number is 
the standard impulse approximation (or extreme closure approximation) which has been widely used throughout the 
literature. The impulse approximation is believed to be a reasonable approximation in intermediate energy nuclear 
physics, i.e. in an energy regime where the kinetic energy of the projectile is large compared to excitation energies of 
the target, however the validity of this assumption has to be tested for individual cases under consideration. 

B. First Order Folding Optical Potential 

In this section we will give the explicit expressions for the 'traditional' first order Watson optical potential. Starting 
from those expressions will lead in a straightforward fashion to an optical potential where the nucleus is treated as a 
composite of clusters. Based on Ec^s. © and (O the first order optical potential as function of external momenta k 
and k' is given by 

(k'|(</.A|PC/P|(/)A)k) = (7eKk',k)= Y. (k'|(</>A|To,;(£:)|0A)k) = (fo.), (12) 

i=N,P 

where £ is the energy of the system. The summation over i indicates that one has to sum over N neutrons and Z 
protons. The structure of Eq. ([T2|) is graphically indicated by Fig. [TJ k^ and k^ are internal variables of the struck 
target nucleon, po and Pg are external target variables. In the following derivation, we will sum over A nucleons for 
clarity of presentation. However it should be emphasized that in practice a sum over N neutrons and Z protons, as 
indicated in Eq. (TTO]) . is done. The energy of the foi operator, £, is a dynamical variable which depends on the total 
energy of the system and the energy of the spectator |I5| similarly to a three-body problem. In this work the common 
approximation to fix it at half the laboratory energy will be used. 

We consider a frame of reference in which kg is the momentum of the projectile, and po is the momentum of the 
target. The individual nucleons inside the target have momenta (ki,k2, ...k^i). Thus 

(kik2k3k4...k^|(/.^) = 5(ki + k2 + kg + k4... + k^ - Po) (C1C2C3C4. .-Ca-I I^Z-a) , (13) 

where the delta function determines the conservation of the absolute momentum for the center of mass (cm.) frame 
of the nucleus, and Q represent the relative momenta of the individual nucleons in the target. The reason for using 
relative momenta is that the wave functions are naturally expressed in such a basis to manifestly express Galilean 
invariance, so any input into the theory will be in terms of the internal momenta ^i and not k. In first order, which is 
the concern of this work, we only need the momentum of the struck nucleon, namely ki = Ci + ■§ ■ Changing integration 
variables from absolute to relative momenta, only the single particle density matrices, /3(CijCi)i ^-re employed. They 
are used to describe the dependence of one particle's relative motion to the remaining (A-I)-core, 

»A-1 „A-1 

p(c;,Ci) = / n < / n ^Cj(0AiciC2C3a-c;i-i) (CiC2C3C4...a-ii</'A). (m) 

•^ 1=2 •' j=2 

The NN foi-matrix can always be written in relative coordinates as 

(qoq'i|7'oi(f)|qoqi) =(5(qo + qi -qo-qi)( ^(^o - ^i)\roi{^)\-j^{<iQ - ^x)j , (is) 



where the (5-function indicates the momentum conservation of the two-nucleon pair. This leads to the following 
expression for the optical potential of Eq. ([T^ . 

(roi) ^ J<i J ^Ci ^(k' + p^ - k - po) 



MCi-Ci-^(k-k')), (16) 



A 
which explicitly gives the following relation 

Cl-Ci = ^(k-k'), (17) 

relating the relative variables C directly to the external variables. For convenience, in the practical calculation, the 
definition 

K^i(k' + k) 

q^k'-k = -^(Ci-c;) 

P^^(Ci+Ci), (18) 

where q is the momentum transfer and K is orthogonal to it, is used. After a series of variable transformations, which 
are outlined in Appendix [Xj one obtains for folding single-particle optical potential 

£/.i(q,K) = 

i—n,p ^ ^ / / \ / 

The above expression shows that one has to carry out a three-dimensional integration over the NN t-matrix and the 
nuclear density matrix. We are employing Monte Carlo methods for the actual computation. 

C. Cluster Model for the Target Nucleus 

Halo nuclei exhibit the distinctive feature of a usually tightly bound core and loosely bound valence nucleons. For 
exploring bound state properties a cluster model has been quite successful (see e.g. [l|, [ia|). Here we propose to 
employ a cluster model for the optical potential of ®He, and thus view ®He as a cluster of a tightly bound a-particle 
and two valence neutrons. As explained in the previous sections, the first order folding optical potential uses the single 
nucleon density as a basic ingredient. That means, in the standard formulation of the optical potential, one deals 
with one active nucleon at a time, and then sums over all active nuclei. This paradigm can be naturally extended to a 
cluster model. The only additional consideration which needs to be taken into account is that now there is an intrinsic 
motion between the different pieces of the cluster. In order to accommodate this let us define a Jacobi momentum 
Pj, representing the relative momentum of the active cluster and the spectator clusters, as 

V3^ = ^{As^V^-A,^,,), (20) 

where the index i characterizes a particular cluster and the index Si represents the spectators for the ith cluster. 
Since this is a Jacobi momentum it is invariant in all frames. The underlying belief is that there is a strongly peaked 
active cluster momentum, centered about pi, which is different from the other spectator cluster momentum, Psj. If 
the nucleons are assumed to be largely independent of each other, as in a single particle picture, then pj . — 0, and 
the single-particle optical potential re-emerges. 

With these additional momenta we can define a correlation density similar to the traditional density of Eq. (|14p : 

Pcorr(Pji,p/i) = / X{dpj[ / W dpj^(0A|pj;pj^...pj'^J (pj^PJ2...pj^J(/)yl), (21) 

•^ 1=2 "^ m=2 



where Nc is the number of total clusters in the target. For the ^He system studied here Nc — 3. The particular view in 
this model is that nucleons within a specific cluster move with the same cm. momentum, which in turn is correlated 
with the momentum of the spectator clusters. This product of densities must conserve the overall momentum in the 
intrinsic frame of the nucleus. It further follows that the sum of all relative momenta in the cluster must be zero. In 
this work, we choose to work with the three-body cluster orbital shell model approximation (COSMA) density [17l[l8|. 
In order to match the choice of momenta to those of the optical potential of the previous subsection, we define 



^.. = '' "'\ (22) 



which is similar to the definitions of P and K defined in Eq. (|T8l) . Thus, the new cluster optical model provides an 
additional sum over the number of clusters 



1 AT ,; ^ ;*, "^ 



c—l,Nc i—ric.p. 

*o.(..i(^K-p)..) ,„(p-±^,,P + ±^<,), ,23, 

where each cluster now defines its own optical potential. With this, the optical potential for ^He consists of two pieces 
as indicated in Fig. [2l an optical potential for the a-core and one for each of the two neutrons, both linked by the 
correlation density between the clusters, 

'"°C/eKq, K) = C/„ + 2t/„ = 
Y. j dV dV,^ PcorriPjJ TO. (q, ^ (^K - p) , ^ ) Pa. (p - ^q, P + ^q 

+2jdP dVj^ PcorriVjJ TO. (^q, ^(^^K - P),£^ p„ (p - ^^q, P 



A-1 



(24) 



The cluster optical potential involves now a six-dimensional integration, which we again carry out with Monte Carlo 
methods. In addition, care must be taken to evaluate each cluster in the cm. frame between projectile and cluster, 
since the scattering takes place in intrinsically different frames as dictated by the variable Vj. However, the final 
optical potential for the nucleus ^He must be evaluated in the cm. frame of the target -|-projectile system. Thus, 
the initial step in the calculation involves first the evaluation of each cluster optical potential separately, and then 
boosting each to the target+projectile system. The explicit details of these transformations are given in Appendix IbI 
Furthermore, it is important to note that each cluster optical potential contains the correlation density. Through 
its variable Vj the momenta of the active particles are constrained. An explicit derivation of the correlation density 
for the ^He nucleus is given in Appendix [Cl 

III. RESULTS AND DISCUSSION 

In this section we evaluate the differential cross section and the analyzing power for elastic scattering of "'He and 
^He using non-local optical potentials in first order in the Watson multiple scattering expansion. Specifically, we 
want to test the influence of the cluster formulation presented in Section II. C on those observables compared to 
single-particle optical potentials of the same order in the multiple scattering expansion. We start with considering the 
recent experimental data for ^He at 71 MeV/nucleon [4|, then we continue our investigation at slightly higher energies 
in order to gain some insight on the behavior of the elastic scattering observables as function of projectile energy. 

As a nucleon-nucleon (NN) interaction we use the nonlocal Nijm-I potential of the Nijmegen group [l^, which 
describes the NN data below 350 MeV laboratory energy with a x^ « 1. For the density matrix of ^He the COSMA 
density of Refs. [Hi JJJ is used. This density consists of harmonic oscillator wave functions for the s- and p-shell. The 
parameters are fitted such that ^He has a charge radius of 1.77 fm and a matter radius of 2.57 fm. 

In addition the elastic scattering observables for proton scattering off ^He is calculated. The proton and neutron 
density matrix for ^He is obtained from a microscopic Hartree-Fock-Bogoliubov (HFB) calculation of Ref. [20|, which 
uses the Gogny DIS finite range effective interaction [2l|. The HFB calculation produces a mean-field potential, 
which in turn is used to modify the free NN interaction inside the nucleus [22, [2J| . This modification of the free NN 
t- matrix has proved to be important in the description of closed shell nuclei at projectile energies below 150 MeV. 

Let us first concentrate on the scattering of ^He at 71 MeV/nucleon. In Fig.|3]we show the differential cross section 
(upper panel (a)) and the analyzing power (lower panel (b)) for elastic scattering of ^He at 71 MeV/nucleon. For 



more detail we show the differential cross section in a linear scale in Fig. |4l The solid (black) line represents the 
calculation with a single-particle optical potential as outlined in Section II. B based on the COSMA density. On the 
log-scale of Fig. [3] the differential cross section looks reasonably well described. However, the linear scale of Fig. |4] 
reveals that a for the small angles, the single-particle optical potential over predicts the differential cross section. 
The analyzing power stays positive for all angles, similar to the predictions in Ref. [6|. Our calculation based on the 
cluster model using again the COSMA density is represented by the short-dashed (blue) line. The differential cross 
section is not very sensitive to the explicit cluster calculation with the exception of the forward angles, where the 
cluster calculation gives a lower cross section in better agreement with the data. The analyzing power, however, does 
not show any improvement though the minimum is slightly shifted to smaller angles, it stays positive, whereas the 
data indicate a negative sign. The folding for the optical potentials with the COSMA density are carried out with 
the free NN t-matrix as input. As has been shown in Ref. [2J| for a variety of heavier nuclei, at energies lower than 
^150 MeV projectile energy, the free NN t-matrix experiences a modification due to the nuclear medium, which can 
be treated as an additional force represented by a mean field acting on the two active nucleons during the scattering 
process [2j|- The dash-double-dotted (green) curves in Figs. [3] and |4] include a modification of the free NN t-matrix 
through an HFB mean field for the alpha cluster only. One advantage to this cluster paradigm is that the calculation 
can utilize the infiuence of a mean field on the free NN interaction where it is most appropriate, i.e. for the strongly 
bound a-core. The results of this calculation produces an analyzing power, which turns negative at 60° and captures 
the shape of the last two measured angles. However, it still over predicts the measured analyzing power at the smaller 
angles. The overall shape of the differential cross section is not modified. Only for the very forward angles the cross 
section is slightly lowered compared to the cluster calculation with the COSMA density as can be seen in Fig. 2] 

It is further instructive to investigate the importance of the correlation density pcorriVjJ in Eq. (PT|) for the cluster 
optical potential. This can be done by realizing that setting PcorrO^j ) = 1 omits the correlations. In Figs. [3] and |3] the 
dash-double-dotted (green) curve represents the calculation based on the cluster formulation whereas the short dashed 
(pink) line represent the same calculation with the correlation density set to 1. The effect on the analyzing power, 
Fig. 131 is relatively small. However, the differential cross section for small angles. Fig. 21 shows visible sensitivity. 
Indeed, one can conclude, that the lowering of the differential cross section for the forward angles is dominated by 
the influence of the correlation density. 

The upgrade of the RIKEN facility will in principle allow to measure the angular distribution of the analyzing 
power for the elastic scattering of ^He at somewhat higher energies. Thus we want to investigate the predictions for 
the elastic scattering observables as function of the projectile energy using a cluster ansatz for the optical potential 
for ^He. As test energies we choose 100 MeV and 200 MeV/nucleon. First, we show in Figs. [5] and [S] the angular 
distributions of the differential cross section and the analyzing power for proton scattering off '*He as function of 
the momentum transfer. The calculations for the projectile energy of 200 MeV shows that here an optical potential 
description of the scattering process is quite good up to 2.5 fm^^ for both, differential cross section and analyzing 
power. In Figs. [S] and [5] two calculations are shown: For the solid line a folding calculation of the optical potential 
has been carried out using the free NN t-matrix, whereas for the dashed line a NN t-matrix modified by a HFB 
mean field has been used. At 200 MeV a folding with the free NN t-matrix is adequate. At the two lower energies, 
a single scattering optical potential describes the differential cross section only up to a momentum transfer of about 
1.75 fm~^. For higher momentum transfers multiple scattering contributions are expected to become important. 
However, only for a deuteron target multiple scattering contributions have been investigated systematically [2^ l2q | 
over a wide range of projectile energies. Thus, we can only speculate about the increasing importance of multiple 
scattering contributions at higher momentum transfer. The analyzing power at 71 MeV projectile energy is not well 
described for small momentum transfers. This is very likely the reason for our over-prediction of the analyzing power 
of ^He for small momentum transfer at this energy [27[. 

In Figs. [7] and [5] we present predictions for elastic scattering of ^He at 100 MeV and 200 MeV/nucleon. The solid 
lines represent the calculations with a single-particle optical potential, whereas the short-dashed and dash-double 
dotted lines are based on the cluster formulation discussed in Section II. C. The short-dashed line uses a free NN 
t-matrix for every piece of the optical potential, whereas the dash-double-dotted line incorporates a NN t-matrix 
modified by an HFB mean field potential for the a-core. Similar to the '^He calculations, at 200 MeV/nucleon the 
modification of the NN t-matrix becomes irrelevant. As discussed for the 71 MeV/nucleon calculations, the cluster 
model lowers the predictions for the differential cross section for small momentum transfers (angles). This feature 
seems independent of the employed projectile energy. For the analyzing power at small momentum transfer the 
contribution of the valence neutrons is not large enough to change the contribution of the a-core part of the optical 
potential. For the 100 MeV calculation the modification of the free NN t-matrix through the nuclear medium is 
still quite pronounced at ~2 fm~^ momentum transfer, whereas the 200 MeV/nucleon calculations exhibit very little 
difference between each other. In Fig. [S] we show the spin-rotation function Q for elastic scattering of ^He at the same 
energies. Obviously, this observable has no chance of being measured. However, since it is an independent observable, 
we consider it instructive to study if it exhibits differences similar to the ones in the analyzing power between a 
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cluster paradigm for the optical potential and a single-particle optical potential. It is interesting to note, that Q is 
very insensitive to any of the modifications introduced, even at 71 MeV/nucleon. 

In order to investigate if the difference of the observables for elastic scattering of "^He and ^He changes with increasing 
energy, we plot in Figs. [TUl and [TT] the differential cross sections and analyzing powers for the two nuclei together at 
71 MeV and 200 MeV as function of the momentum transfer. One could speculate that when employing a cluster 
model for the ^He nucleus, at higher energies the piece of the optical potential due to the a-core might dominate the 
observables (possibly the analyzing power). However, this does not seem to be the case, the behavior of the differential 
cross section is very similar at the two energies. Even the shift of the minima in the analyzing power between the 
two nuclei is the same for both energies, namely roughly 0.25 fni^^. Obviously, experimental information will have 
to decide if a cluster ansatz as presented here captures a major part of the underlying physics. 



IV. SUMMARY AND OUTLOOK 

In this work we introduced a cluster formulation for the optical potential for calculating scattering observables for 
elastic scattering of three-body halo nuclei. We first reviewed a traditional single-particle full-folding optical potential 
in first order in the Watson multiple scattering expansion, and then showed how one can naturally extend this optical 
potential to introduce the cluster structure of a halo nucleus. Here we concentrate on the ^He nucleus. However, the 
formulation we introduced can be further extended to four or five-body clusters, e.g. to ^He. For our calculations 
we used the density matrix of the three-body cluster orbital shell model approximation (COSMA) introduced in 
Refs. [13, ll8i for the ^He nucleus. This density matrix is based on single harmonic oscillator wave functions for the s- 
and p-shell of ^He and allows a straightforward calculation of the required correlation densities needed for the optical 
potential. The resulting folding optical potential contains a six dimensional integration over internal vector momenta, 
which is calculated via Monte Carlo integration. 

We calculated the angular distribution of the differential cross section and the analyzing power at 71 MeV, 100 MeV, 
and 200 MeV/nucleon and compared our results with experimental data at 71 MeV/nucleon [J, [2^ . We find that the 
cluster model lowers the cross section for the small angles and brings it closer to the data. Though we do not describe 
the very small analyzing power at the small angles, we find that the cluster formulation together with a 'hybrid' 
ansatz, in which the optical potential for the alpha-core is calculated with a NN t-matrix modified by a mean field of 
the alpha is able to produce a negative analyzing power at larger angles as suggested by the data. Our predictions for 
the higher energies indicate that the lowering of the differential cross section for small momentum transfers (angles) 
using a cluster paradigm remains visible. The cluster ansatz for the optical potential continues to predict a negative 
analyzing power at larger momentum transfer. Eventually experimental information should be become available to 
see if these predictions capture the bulk of the physics of the reaction at higher energies, or if there are additional 
theoretical pieces necessary to understand this reaction, preferably as function of scattering energy. 



Appendix A: First Order Full-Folding Optical Potential 

In this appendix we will give explicit steps for arriving at the 'traditional' first order Watson optical potential of 
our calculations. The complete derivation is given in Ref. [15'], however for the convenience of the reader we want to 
give a shorter summary here. 

Using the (5- function of Eq. (J13p . which determines the conservation of the absolute momentum for the center of 
mass (cm.) frame of the nucleus and the relative momenta of the individual nucleons in the target, we obtain for the 
optical potential of Eq. (IT^ 

(foi) = (k'|(0^|fo,(£:)|0A)|k) 

= /n^^- IW^i (0AiCiC^aci-c;i-i)'5(p'-p;,) (k'kiifoi(£)ikki) 
•' j=i ■^ 1=1 

A 

l[5m~k,)5{p-po) (C1C2C3C4...CA-1I0A), (Al) 

i=2 

where p = J2i=i ^i, P' = Si=i ^i- Without losing generality, we consider for now nucleon '1' being the active nucleon 
in the target. The additional delta functions arises because the target nucleons are independent of tqi. A Galilean 



invariant choice of internal variables are the Jacobi coordinates, 

k.^a + f 

, > Ci , p 

k2 = C2 - 



ks - C: 
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A-l A 
Ci _ C2 P 

A-l A-2 A 



In first order, which is the concern of this work, we only need the momentum of the struck nucleon, namely ki = Ci + 5 • 
Changing integration variables from absolute to relative momenta in Eq. (JA1[) results in 

.A-l ^A-1 



roi) - / n '^cjrfp' / n ^cirfp ('/'AiciC2C3c;-CA-i)'^(p' - Po) 

^ A 1 ^ 1 1 



j=l •' 1=1 

(k',Ci + 4|roi(f)|k,Ci + -^) n ^(Cj - Cj)-? ( ^p' -Ci- ^P + Ci 



A-l 

n 

J=2 



<5(p-Po)(CiC2C3C4...Ca-i|<^a). (A3) 

Eq. (|A3[) indicates that one only needs the single particle density matrices p(Ci, Ci)> describing the dependence of one 
particle's relative motion to the remaining (A-l)-core, and given in Eq. ([T4)). Inserting p(Ci,Ci) into Eq. (|A3[) and 
evaluating S{p — po) and S{p' — Pq) leads to 



ifoi) 



JdC JdCik' Ci + ^ |roi(f)|k Ci + ^) p(Ci,Ci) 



^Po-Ci-^Po + Ci). (A4) 



Inserting the NN foi-matrix of Eq. (eq:2.8.1c) and taking advantage of its conservation of the cm. momentum 
Eq. (jA4p becomes 

(roi) = J<i I dCi S{-k' + p^ - k - Po) 

(^(k'-Ci-^)|roi(f)|i(k-Ci"^)) piCiXi) 

'5(^p[,-Cl-^Po + Ci). (A5) 

The first delta function describes the overall momentum conservation, which can be used to reduce the integral of 
Eq. (jASP to three dimensions, i.e. the integral given in Eq. (|T6l) . In our practical calculations we use the variables q, 
K, and P, which are given in Eq. ([T8| . The inverse relations are given by 

k = K iq 
k' = K+iq 

Cl=P-^q. (A6) 

Substituting those variables into Eq. (IT5|) leads to 
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Here we dropped the overall momentum conserving delta function, which is carried out when evaluating the cross 
section. Since the preceding derivation based on the projectile particle '0' and active target particle '1' is general, 
and thus can be repeated for all N target neutrons and Z target protons, one obtains as final expression for the 
'full-folding' optical potential 

C/e,(^K+lq,K-iq^ = 

f^ A-1 A-1 \ , . , 

P.(P-^^q,P + ^^qj. (A8) 

From Eq. (jA8[) we can read off the momenta of the NN t-matrix as 



ifK-P + ^q-4 
2 V 2A ^ A 



kW = ;^(K-P + ^^l^q-^). (A9) 



For numerical calculations we prefer 



(Inn — ^NN ~ ^NN — k — k — q 
K^^ = I (kW + kA.^) ^ i (^K - P - ^^^) . (AlO) 

Since q is a momentum transfer, it is invariant under Galilean transformations, i.e. q^NN = Q- Rewriting the optical 
potential in terms of q and K gives 

I'/U. (,, i (k - P - £^) ,.) . (P - ^,.P . ±^,) . ,An, 

Rewriting Eq. (JA11[) in variables of the NA cm. system requires that k + po = k' + pg = 0. Thus, the dependence 
of Po drops out, and one obtains for the 'full-folding' single-particle optical potential of Eq. ([T9l) 



Cf.i(q,K) = 

Appendix B: Transforming the optical potential 

We first want to calculate the optical potential for a nucleon scattering off a cluster 'i' in the cm. frame of the 
projectile and the cluster 'i' (e.g. cluster 'i' could be the a particle within the ^He nucleus). Let us use the subscript 
Ci to denote that frame, whereas quantities without subscripts shall be interpreted as given in the A -I- 1 cm. frame. 
The overall conservation of momentum for the cm. cluster frame assumes 



^c. 



kic. ^k'c.+k'ic, =0, (Bl) 



where kc^ is the projectile and ki^;. is a typical target nucleon inside the nucleus, within the ith cluster. 
Starting with the definition of the density given in Eq. (|A2[) in the A 4- 1 frame, 

Ci = ki--, 

this becomes in a specific cluster frame 

Ci=kic, -P^. (B2) 
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The momentum Ci does not carry a cluster subscript since it is defined in the traditional intrinsic frame of the single 
particle density. Eq. (|B2p defines how this intrinsic variable is related via a Galilean transformation to the cluster 
frame. Using the notation of Eq. ([20| . this can be broken up into the active particle and the spectator 

Ci = kic^ -J . (B3) 

In this center of mass frame one has kp. = —pic , where k^; is the momentum of the projectile in the cluster frame. 
Then using Eq. (jB3|) . this result together with understanding that the spectator momentum does not change during 
the collision, we can calculate the difference 

Ci - Ci' = (kic. - k;^,) - i ^'^'~^^n = ^q, (B4) 

remembering that q, the relative momentum transfer, is invariant in all frames. If we allow the same definition for 
the average momentum of the target nuclcon, Eq. (fT8)) 

(B5) 



then the inverse equations of Eq. (|A6 



p 


^^(Ci' + Ci), 


Ci = 


--^o 


Ci' 


-^-'-\ 



(B6) 

have the same form. This makes sense, since q is invariant and P is written in the intrinsic frame of the nucleus. 

The momentum arguments of the foi-matrix are defined in Eq. (jA9p for the yl + 1 frame. However, we also need to 
redefine them in the cluster frame Ci. Rewriting Eq. (|A9|) and using Eqs. (IB3p and (IB6p . we can write the relative 
momentum between the projectile, k^ and the struck nucleon, ki^-. as 

k(7. — Psic 

kc, - kic, = kNNc, = kci ^ Ci H ^ -■ (B7) 



The primed momentum is similarly given as 

k'c^-k'=k'^=k'c,-C^' + ^i-^. (B8) 



The difference between Eq. (|B7|) and Eq. (|B8p gives the momentum transfer q as expected. We can define the sum of 
the momenta of the foi-matrix as the sum of Eqs. (IB7p and (|B8I) as 



Knnc, = ;; = ^n-C^Ci + kcj - 7;(Ci + Ci ) ;;— , (B9) 



remembering that the spectator momentum does not change during the collision. In the spirit of Eq. (J18p we can 
define the sum of the cluster momenta as Kc^ = ^ (kNN Ci + k^Nc-) ^^^'^ then rewrite Eq. (jB9l) as 

Knnc, = 2 ^ A • A~' ^ ' 

Examining Eqs. (|A9p through (|19p. we see that the momentum argument Knn of the fgi-matrix has picked up an 
extra term involving the spectator momentum in the cluster frame. The rational for this is simple, the struck nucleon 
is acted upon in a specific cluster, but the total density contains both cluster and spectators. 

We can rewrite Knn in terms of the correlation momentum Vj ■ by multiplying Knnc- by the number of nuclcons 
in the «th cluster, Ai, 

MK^Mc, = A, Ai^Kc. - AP - A^. (Bll) 

After some manipulation as well as using the definition of Vj ■ found in Eq. (j2ip , one finds 

MKnmc. = {A, + l)Kc, - A,P + Vj^, (B12) 
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so a cleaner definition of the average momentum can be written as 



ifjvivc. = ^1— Kc, -P + ^. (B13) 

Thus, for a specific cluster frame we can write the optical potential as 

In the case of ^He this is the optical potential for a proton on the alpha core (or the proton projectile on one of 
the neutrons). We have not worried about the correlation density transformation since this is based on a relative 
momentum Vj^ and is thus invariant during a frame transformation, as are P and q, the variables of the traditional 
single particle density. 

In order to be useful, Eq. (|B14p must be transformed from each individual cluster frame back to the nucleon-nucleus 
frame, so that it can be summed with the other clusters which make up the target nucleus. Scattering observables 
can then be calculated in the cm. frame of the A + 1 system following Eq. ^1% . As it stands, each cluster has its 
own unique cm. frame, and thus they cannot be summed until they are transformed back to the unique ^ + 1 cm. 
frame. The only argument of concern, because it is not invariant, is the momentum Kp. in the fgi operator. 

Employing conservation of momentum in both frames, we can define how the cluster frame relates to the nucleon- 
nucleus ^ + 1 frame. Setting relative velocities equivalent in the two different frames leads to 

A,k - Pi = v4,kc. - Pic. = {A^ + l)kc. , (B15) 

where the last equivalence is given because in the cm. of the cluster frame kc^ = — pi^.. Another second relation 
between the two frames is gained by examining the Jacobi momentum of Eq. ((20)) in the A + 1 frame, 

PJ^ = ]4(^«P»-^«Ps^)• (B16) 



Rearranging gives 



^ (p. + PsJ + P,,, = -^k + p,., (B17) 



where again the last equivalence is found by recognizing that in the A + 1 frame the cm. momentum is defined as 
k + Pi + Psj = 0. Inserting the result for pi from Eq. (jB17|) into Eq. (|B15P we can after some manipulation compare 
the momenta between the two frames: 

The same relationship can be developed for the primed momentum 

A 



k' 



A- 



-i-((A. + l)k'^^+p;,). (B19) 



Adding these two equations gives 



^ = -A~A ((^* + ^^^^- + ^'^ ■ ^^^°) 



This relation is the transformation prescription for K from the cluster frame to the A + 1 frame. Solving for Kc^ in 
Eq. (jB20p and plugging it into the optical potential in the cluster frame of Eq. (|B14[) we are then able to express this 
potential completely using invariants or nucleon-nucleus A H- 1 variables. 






This is Eq. (12^ from Section II. C for the cluster optical potential in the nucleon-nucleus frame. 
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Appendix C: Correlation Density for the Cluster Approach 

The cluster approach developed in this work uses a correlation density relating the clusters, which is given in 
Eq. dH]), 



"(Pjl'Pji) =y IlMlJ n ^Pjm(<^^IPjiPj2-PjArj (PJlPJ2-PJArJ'^^)' (d) 



No „ N. 

/ t \ I 

1=2 '' m=2 

This density correlates the momenta between the various clusters and elevates the this approach beyond the indepen- 
dent single particle picture. For the explicit derivation, let us start from definition given in Eq. (|20|) 

Pj»^^(^«P*-^^Ps.)- (C2) 

This definition is invariant from the frame of consideration. Thus it can be applied in the frame of the intrinsic 
density, where it describes the difference in momenta between cluster i and its analogous spectator particles, Psj, in 
the laboratory frame. In this same frame the total momentum between active and spectator particles should add up 
to zero, at least before the collision. Thus, Pi = — Ps^, and therefore Pj. — Pi and Pj/ = pi — q. Again, in the 
intrinsic frame of the nucleus one has 

Vn = ^^^ = P.,: - \ (C3) 

For the ^He nucleus consisting of three clusters, we can write the correlation density as 

rr(Pji,Pji) = 

d3p,^^d(z-pi) d^p^^^ d3p,^2 d3p;^2 *(P«,l) *(P.,2) *(P^a) *(Pi.2) 
/corr(f^Si,f^S2)(5(Pi -Psii -iPs,2) '^(Psil -Psa) '^(P«i2 ~Psi2)' (^4) 

where the two spectator momenta are labeled '^sn and Psi2- The integration variables are over the momenta (before 
and after the scattering) of the two spectators. We also integrate over the relationship of the active particle's momen- 
tum to the quantization axis of ^He. The first momentum conserving delta function preserves the cm. momentum 
(the last spectator therefore needs not be integrated over), the momenta of active and spectator particles must add 
up to zero. The remaining two delta functions require that the momentum of the spectator clusters do not change. 
In practice, Eq. (jC4p is for ^He a four dimensional integration. In this work, the wave functions are the single particle 
wave functions, where the momentum is converted to single particle form by simply dividing by the cluster mass. In 
one had a density defined using both in a cluster and a single particle paradigm, then this formulation would allow 
for increased dynamical detail. The angular correlation function, fcom gives the correlation weighting assuming that 
the spectator neutrons are in a given orbital shell, where (llsj,!!^^) are the solid angles subtended by the spectator 
valence nucleons, in this case the p| shell. 

For calculation the correlation density we assume, that the two valence neutrons are in the p| shell, and the total 
angular wave function can be written as 



/■ 



\{i-Pi2)yp{<:i)y! MC2) + ^( 



^ps = 7t(i - Pr2)ypi':i)y^ MC2) + ^^(i - Pi2)ypiCi)y? MC2), (cs) 



where y^ ' (C) are the traditional spin spherical harmonics and C, is the single particle intrinsic momentum. The 
antisymmetric of the wave function with respect to the two neutrons is given by the operator (1 — P12). The correlation 
function, fcorr{^sn ^s^), defined in Eq. ()C4p can be with the help of Eq. (jCSP defined as 

fcorr = %A(lX2)^pA^lX2)- (C6) 

2 ^ 

This defines the angular probability for the two spectator neutrons when the alpha core is active. This local correlation 
function is also used as the approximate probability when a neutron is the active particle, i.e. when the full calculation 
is in fact off-shell. The alpha core has an unweighted angular distribution, approximated by the COSMA density, as 
being completely in the s-orbital. 

Explicitly inserting the spherical harmonics into Eq. (jC5l) . we obtain 

1 2 

Icorr = -5:5-^ (2 COS 6*1 cos 6*2 - sin 6*1 sin 6*2 cos((/)i - ^2)) 



327r2 



327r 



( sin^ 6*1 cos^ 6*2 + cos^ 9x sin^ 62 - sin 26*1 sin 26*2 cos(0i - 02 )) ■ (C7) 
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The first term is the result if the total spin projection of the two neutrons added up to zero and the second term is 
due to the total spin projection being one. 

Once the four dimensional integral for the cluster correlation, Eq. (jC4|) . is calculated, it is normalized to one and 
then used to augment the definition of the optical potential. De facto, it constrains the momentum of the cm. of the 
active cluster relative to the spectators. 
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FIG. 1. Diagram for the standard optical potential matrix element for the single-scattering approximation. 
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Pi p^c-fo. 




FIG. 2. Diagram for the cluster optical potential for ^He based on the single-scattering approximation. 
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FIG. 3. (color online) The angular distribution of the differential cross section,^, (upper panel (a)) and the analyzing power, 
Ay, (lower panel (b)) for elastic scattering ^He at projectile energy 71 MeV/nucleon as function of the cm. angle. The 
calculations are performed with optical potential obtained from the Nijmegen I potential [13| for the NN interaction. All 
optical potential are folding, non-local optical potentials described in the text. The solid line (black) represents the calculation 
based on a single-particle optical potential employing the COSMA density of Ref. [ITI ]. For the short-dashed line (blue) the 
cluster ansatz together with the COSMA density is used. The dash-double-dot line (green) represents a calculation based on 
the cluster formulation, however the NN t-matrix for the core optical potential is modified by a mean field obtained from a 
HFB [20l . [2H calculation. The short dashed line (pink) represents the same calculation, but neglects correlation of the clusters. 
The data are taken from Refs. P. |28J. 
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FIG. 4. (color online) Same as Fig. 3 except that the angular distribution of the differential cross section (^) is plotted with 
a linear scale. The data are from Ref. 1281. 
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FIG. 5. (color online) The angular distribution of the differential cross section (^) for elastic proton scattering off He at 
projectile energies 71 MeV, 100 MeV, and 200 MeV as function of the momentum transfer. The calculations are performed 
with optical potential based on the Nijmegen I potential [IQI for the NN interaction. For the a-core an HFB density according 
to Refs. 20| . [2ll | is employed. The solid lines (black) show the calculations based on the free NN t-matrix, while the dotted 
(red) lines are based on calculations modifying the NN t-matrix with a mean-field consistent with the HFB Q-core. The data 
are from Refs. i29- 32!l . 
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FIG. 6. (color online) The angular distribution of the analyzing power for elastic proton scattering off *He at projectile energies 
71 MeV, 100 MeV, and 200 MeV as function of the momentum transfer. The meaning of the curves is the same as in Fig. 5. 
The data are from Refs. j29-32]. 
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FIG. 7. (color online) The angular distribution of the differential cross section (^) for elastic scattering of ®He at projectile 
energies 71 MeV, 100 MeV, and 200 MeV/nucleon as function of the momentum transfer. The calculations are performed with 
optical potential based on the Nijmegen I potential [IQI for the NN interaction. For the solid (black) line the COSMA [131 has 
been used as single-particle density. The short-dashed (blue) line incorporated the cluster structure into the optical potential 
using the COSMA density for all clusters. For the dash-double-dot (green) line the free NN t-matrix has been modified with 
the HFB mean field. The data are taken from Ref. [1. [28l|. 
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FIG. 8. (color online) The angular distribution of the the analyzing power (Ay) for elastic scattering ^He at projectile energies 
71 MeV, 100 MeV, and 200 MeV/nucleon as function of the momentum transfer. The meaning of the lines is the same as in 
Fig. 7. The data are taken from Ref. [4, 28]. 
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FIG. 9. (color online) The angular distribution of the the spin rotation function Q for elastic scattering ^He at projectile 
energies 71 MeV, 100 MeV, and 200 MeV/nucleon as function of the momentum transfer. The meaning of the lines is the same 
as in Fig. 7. 



22 



10^ b- 

a 

"^ 2 



10" 



10^ b- 

^10° 

a 

|io-^ 

10^ 





1 ' 1 ' 1 ' 1 ' 1 E 


1 lllllll 1 lllllll 1 lllllll 1 lllllll 1 Illll 


'^•.^***'». ("He 71 MeV) n 
■^ . . " " ■ - -"1 

(^HeTlMeVr*^ "^ 

1 1 1 1 1 1 1 1 1 ~ 



10" 



: » 1 
= \ 


1 


1 

1 


*-= 


1 

1 


1 

1 


1 ' 1 ' 1 1 

r^"-... (^HeZOOMeV) 1 
"*^. ' " * • . 

(^He200MeV) '^••i 

1 1 1 1 1 " 


1 1 1 1 1 III 


1 





0.5 






1 




1.5 2 2.5 



q (fni ) 



FIG. 10. (color online) The angular distribution of the differential cross section, ^, for elastic proton scattering off *He 
(short-dashed red) and ®He (dash-double-dot green) at projectile energies 71 MeV (upper panel (a)) and 200 MeV/nucleon 
(lower panel (b)) as function of the momentum transfer. Both calculations are performed with optical potential based on the 
Nijmegen I potential [ly] for the NN interaction, which for '*He (in the ^He case the *He-core) is modified by the HFB mean 
field. For the ^He calculations the single particle nucleons are described by the COSMA density, and the cluster ansatz is used. 
For the *He calculations the HFB density ,2^, i21,] is employed. 
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FIG. 11. (color online) The angular distribution of the analyzing power for elastic proton scattering off ^He (short-dashed red) 
and ^He (dash-double-dot green) at projectile energies 71 MeV (upper panel (a)) and 200 MeV/nucleon (lower panel (b)) as 
function of the momentum transfer. The meaning of the curves is the same as in Fig. 1101 



